Monte Carlo simulations of the four-dimensional XY spin glass at low temperatures 
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We report on the results for simulations of the four-dimensional XY spin glass using the parallel 
tempering Monte Carlo method at low temperatures for moderate sizes. Our results are quali- 
tatively consistent with earlier work on the three-dimensional gauge glass as well as three- and 
four-dimensional Edwards- Anderson Ising spin glass. An extrapolation of our results would indicate 
that large-scale excitations cost only a finite amount of energy in the thermodynamic limit. The 
surface of these excitations may be fractal, although we cannot rule out a scenario compatible with 
replica symmetry breaking in which the surface of low-energy large-scale excitations is space filling. 
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PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

There has been an ongoing controversy regarding the 
spin-glass phase. There are two main theories: the 
"droplet picture" (DP) by Fisher and Huseo and 
replica symmetry breaking picture (RSB) by Parisi 
While RSB follows the exact solution of the Sherrington- 
Kirkpatrick model and predicts that excitations which 
involve a finite fraction of the spins cost a finite energy 
in the thermodynamic limit, the droplet picture states 
that a cluster of spins of size I costs an energy propor- 
tional to l e , where 9 is positive. It follows that in the 
thermodynamic limit, excitations that flip a finite clus- 
ter of spins cost an infinite energy. In addition, the DP 
states that these excitations are fractal with a fractal di- 
mension d s < d, where d is the space dimoasion, whereas 
in RSB these excitations are space filling,Q i.e., d s = d 

Krzakala and Martinja as well as Palassini and Youn 
(referred to as KMPY) found, on the basis of numerical 
results on small systems with Ising symmetry, that an 
intermediate picture may be present: while the surface 
of large-scale excitations appears to be fractal, only a 
finite amount of energy is needed to excite them in the 
thermodynamic limit. In the context of their work, it is 
necessary to introduce two exponents, 9 and 9', where L 
is the typical energy for an excitation induced by a change 
in boundary conditions in a system of linear size L, and 
L describes the energy of thermally excited system-size 
clusters. Subsequently, similar xesults were found for the 
three-dimensional gauge glass,Q which has a continuous 
symmetry but is known to have a finite T c . 

The differences between DP and rRSB|-can be quan- 
tified by studying the distributiona'clofiyO P(q) of the 
spin overlap q defined in Eq. (Q) below. For finite sys- 
tems, the DP predicts two peaks at ±<7ea, where <?ea is 
the Edwards- Anderson order parameter, as well as a tail 
down-to q = that vanishes in the thermodynamic limit 
likeBB ~ L~ 9 . On the contrary, RSB predicts a non- 
trivial distribution with a finite weight in the tail down 
to q = 0, independent of system size. 

Earlier work that studied the nature of tiie-jSpin-glass 
state has focused on the Ising spin glass JatffirtJ though 
some work has also been carried out on the gauge glass 



model of the vortex glass transition in superconductors.0 
Here, we consider a vector spin-glass model, the four- 
dimensional XY spin glass, .-which is known to have a 
finite transition temperatureOJ T c with T c ~ 0.95. We 
perform Monte Carlo simulations for a modest range of 
sizes down to low temperatures. IT ~ 0.2T C ) using the 
parallel tempering Monte CarloE_3li3 technique. Our main 
result is that that P(0) does not appear to decrease with 
increasing system size for the range of sizes studied. 

We also look for information on the surface of the large- 
scale low-energy excitations by studying the "link over- 
lap" defined in Eq. ( |l3| ) below. The data for this quan- 
tity suggests that the surface may be space filling, i.e., 
d s = d, as in RSB, though the small range of sizes pre- 
cludes us from making a firm statement on this and a 
scenario compatible with the DP is also viable in which 
d s < d. 

The layout of the paper is as follows: In Sec. II we 
describe the model and the measured observables. We 
discuss our equilibration tests for the parallel tempe ring 



Monte Carlo method for this specific model in Sec. |H 
Our results are discussed in Sec. IV. Section |y| summa- 



rizes our conclusions and presents ideas for future work. 



II. MODEL AND OBSERVABLES 

The XY spin glass consists of two-component spins of 
unit length on a hypercubic lattice in four dimensions 
with periodic boundary conditions. The Hamiltonian is 
given by 



H = - Jk 

{hi) 



Si • S 



(1) 



where the sum is over nearest neighbors, the linear size 
is L, the number of spins is N — L 4 , and Si = (Sf, Sf) 
is an XY spin. Since |Sj| = 1, one can parametrize the 
spins as = [cos(<^j),sin(</>j)] with fa G [0, 27r]. The 
Hamiltonian then transforms to 



H 



5> 



(2) 
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The Jij are chosen according to a Gaussian distribu- 



tion with zero mean and standard deviation J, i.e., 



2nJ 



exp 



7 2 
2J 2 



(3) 



Unless otherwise stated we will take J = 1. 

We concentrate on two observables, the spin overlap 
q and the link overlap qi. The (tensor) spin overlap is 
defined in terms of the spin configurations of two copies 
of the system, denoted by (1) and (2), as follows: 



1 N 



(l)o(2) 



,M, v G {x,y}. 



(4) 



In analytic work, the spin-glass order parameter is de- 
fined to be the average of the trace of q^ v . To be precise, 
for L — > oo, the order parameter tensor is predicted to 
be of the form 



q/2 






q/2 



(5) 



However, this implicitly assumes that the symmetry has 
been broken by a small field, which is inconvenient to im- 
plement in numerics, so we adopt the following equivalent 
procedure. We apply all possible symmetries (rotations 
and reflection) to one replica and take the largest value of 
the resulting trace. Consider first rotations under which 
/ where 



<i 



Qxx 
Qyx 



Qxy 



(6) 



Maximizing Tr(g') with respect to the relative rotation 
angle between the replicas gives q\ , where 



5i 



(Qxx + qyy) 2 + (Qyx Qxy) z 



(7) 



The rotation also makes the two off-diagonal pieces equal, 

i- e, 5 Qxy Qyx ■ 

We also must consider how the q^ transform under 
reflections of the angles of the spins in one replica, fa — > 
—fa. It is easy to see that under this transformation 
Qi — ► 92 and vice-versa, where 



Q2 = \J (Qxx - Qyy? + (Qxy + Qyx) 2 
= y (qxx - Qyy) 2 + (qxy + %r) 2 , 



(8) 



where the second line follows after some algebra. Since 
the spin-glass order parameter is obtained by maximizing 
the trace with respect to all symmetry transformation, it 
is given by 



q = max{qx,q 2 }. 



(9) 



We use the notation q, somewhat inconsistently, for the 
spin-glass order parameter to conform with notation in 



other work. The spin-glass order parameter function in 
RSB theory, P(q), is given by the distribution of q in 
Eq. (|). 

We also define the smaller of q\ and q 2 by q, i.e., 



q = mia{g 1 ,g 2 }. 



(10) 



If the order parameter tensor tends to the form in Eq. (|5|) 
for L — > oo then q — > in this limit. We shall see that 
our results support this. 

If we are willing to assume that the form in Eq. (|5|) 
applies in the thermodynamic limit then we can obtain 
the-spin glass order parameter distribution a little more 
simply from the quantity 



Q = \Jllx + q 2 yy + qjx + <&yi ( U ) 

which is invariant under symmetry transformations. 
Since 



2Q 2 = q 2 



(12) 



then, if g — ^ for L — > oo, the distributions of q and 
V2Q are the same in this limit. 

The link overlap is defined, quite simply, by 



(13) 



where Nf, — Nd is the number of bonds (d = 4 is the 
space dimension). Since this is already invariant under 
global symmetry operations we do not need to consider 
the effects of rotations and reflections as we did for the 
spin overlap. The link overlap can be expressed in terms 
of spin angles by 



1 



»?») 



(14) 



While a change in q induced by flipping a cluster 
of spins is proportional to the volume of the cluster, 
qi changes by an amount proportional to the surface 
of the cluster. The weight in P(q) for small q varies 
as L~ , where 8' was introduced in Sec. I. In addi- 
tion, we expect the variance of the link overlap to fit 
to a form Var(q/) ~ where, as shown in Ref. ||, 

m = ff + 2(d-d s ). 



III. EQUILIBRATION 

For the simulations-, we use the parallel tempering 
Monte Carlo mcthodJl3'E3 In this technique, one sim- 
ulates identical replicas of the system at Nt different 
temperatures, and, in addition to the usual local moves, 
one performs global moves where the temperatures of 
two replicas (with adjacent temperatures) are exchanged. 
This allows us to study larger systems at lower temper- 
atures than with the conventional Monte Carlo method. 
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TABLE I: Parameters of the simulation. 



iVsamp is the number 
is the total 



of samples, i.e., sets of disorder realizations, 7V S , 
number of sweeps simulated for each of the 2Nt replicas for 
a single sample, and Nt is the number of temperatures used 
in the parallel tempering method. 



L 


■Afj s amp 


iV sweep 


N T 


3 


1 x 10 4 


3.0 x 10 4 


39 


4 


2 x 10 3 


4.0 x 10 5 


39 


5 


1 x 10 3 


2.0 x 10 6 


39 



Since we require two copies at each temperature to de- 
termine the spin and link overlaps, see Eqs. (^) and (13), 
we actually simulate 2Nt replicas. 

The lowest temperature has to be far below T c ~ 0.95 
and yet high enough that a range of sizes can be simu- 
lated. We chose the value of 0.2. The highest temper- 
ature has to be such that the system equilibrates very 
fast, and we chose 1.498. The intermediate temperatures 
arc determined empirically provided that the acceptance 
ratios of the moves interchanging the replicas are larger 
than about 0.4 and are all roughly equal. 

Table | lists the parameters of the simulation; N samp 
(number of samples), iV sweep (total number of sweeps per- 
formed by each set of spins), and Nt (number of tem- 
perature values). 

It is important to ensure that the system is equili- 



brated. However_the equilibration test proposed by 
Bhatt and YoungcJ does not work with parallel tem- 
pering Monte Carlo because the temperature of each 
replica does not stay constant throughout the simula- 
tion. .-Here we use the method introduced by Katzgraber 
et alB for short-range spin glasses with a Gaussian distri- 
bution of exchange interactions that relates the average 
energy to the link overlap. By performing an integra- 
tion by parts with respect to Jy of the average energy 
U = [<«)] av (< 0), we obtain 



[<«>]« 



Is 



2T\U\ 



(15) 



where z is the number of nearest neighbors, (■ • •) denotes 
a thermal average, and [• ■ -] av denotes an average over 
the disorder. The quantity q s is given by 



(16) 



where the sum is over pairs of neighboring spins. The 
simulation is started with randomly chosen spins so that 
all replicas are uncorrelated. This will have the effect 
that both sides of Eq. (|lj) are approached from opposite 
directions. Once they agree, the system is in equilibrium 
as can be seen in Fig. |l| for T = 0.2 (to be compared with 
T c ps 0.95),Eil the lowest temperature simulated, and for 
L = 3. We show data for the smallest size since it allows 
us to generate more samples for longer equilibration times 
to better illustrate the method. For larger system sizes 




sweep 



FIG. 1: A plot of [(qi)] av (the link overlap), qi(U,q B ) de- 
fined to be the RHS of Eq. (^), [(?)]av the spin overlap, 
and [(g)] av defined in Eq. (TTu) , as a function of Monte Carlo 
sweeps iVswccp for each replica, averaged over the last half 
of the sweeps. For equilibration, [(<?i)]av and qi(U,q s ) should 
agree. The two sets of data approach each other from oppo- 
site directions and, once converged, do not seem to change 
at longer times, indicating that the system is equilibrated. 
The data for [(g)]av and [(g)]av show that they too have equi- 
librated in roughly the same equilibration time. While not 
shown here, data for higher moments of the different observ- 
ables have the same equilibration time as the link overlap 
[{<7i)]av. (Data for L = 3, T = 0.2, and 3230 samples). 



we stop the simulation, once the data for [(qi)] av an d the 
right-hand side (RHS) of Eq. ( |l~5| ) agree. 

Because the XY spin glass has a vector order param- 
eter symmetry, we discretize the angles of the spins to 
N^ — 512 to speed up the simulation. This number is 
large enough to avoid any crossover effects to other mod- 
els as discussed by Cieplak et al&3. To ensure a reason- 
able acceptance ratio for single-spin Monte Carlo moves, 
we choose the proposed new angle for a spin within an ac- 
ceptance window about the current angle, where the size 
of the window is proportional to the temperature T. By 
tuning a numerical prefactor, we ensure the acceptance 
ratios for these local moves are not smaller than 0.4 for 
each system size at the lowest temperature simulated. 



IV. RESULTS 

Figures | and | show data for P(q) for T = 0.20 and 
0.42, respectively. In both cases we see a peak at large 
q and a tail for smaller q that does not extend to q = 0. 
However, it is not surprising that there is a "hole" at 
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FIG. 2: Data for the spin overlap distribution P(q) at tem- 
perature T — 0.20 for different system sizes. Note the loga- 
rithmic vertical scale. The lines go through all the data points 
but, for clarity, only some of the data points are shown. 
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FIG. 3: Same as Fig. | but at temperature T = 0.42. 

small q since q is defined to be the maximum of q\ and 
52- If q = min{qi,(72} tends to zero at large L, which 
is expected as discussed above, then, in RSB theory, the 
tail would extend to smaller values of q for larger L while 
maintaining the same height. Looking at Figs. || and [|, 
this seems to be the case, at least for the range of sizes 



FIG. 4: Data for the overlap distribution P(q) at tempera- 
ture T — 0.20 for different system sizes. The weight in the 
distribution tends towards q = for increasing L. 



that we have been able to study. 

In Fig. | we show data for P(q) at T = 0.20. As 
expected, the distributions seem to collapse to zero for 
increasing system size. Figure |^ shows the variation of 
the mean of q with L on a log- log plot. The data have 
been fitted to-straight lines with slopes shown. The qual- 
ity of the fitsta is only moderate; Q = 0.06, 0.09, and 0.04 
for T = 0.200,0.247, and 0.305, respectively. Given the 
rather small range of sizes, and hence the likelihood of 
systematic corrections to scaling, we feel that the data 
are consistent with [(?)] a v - * for L — > oo. Since q > 0, 
if [(g)]av = then the whole distribution collapses to 
q = 0. 

Lastly we present in Figs. || and [?] results for the distri- 
bution of the link overlap P{qi). There is a pronounced 
peak at large qi values as well as the hint of a shoulder 
for smaller values in the T = 0.20 data. The width of the 
distribution decreases with increasing system size. This 
is demonstrated in Fig. ||[ which shows the variance of qi 
against system size L for several low temperatures. 

There is some curvature in the data for Var (<?/), so first 
we attempt a three-parameter fit of the form 

Var(qi) = a + bL~ c , (17) 

finding small but finite values for a, see 

Table ||. As we 

have the same number of data points as variables, we 
cannot assign fitting probabilities to the fits. We also 
attempt a power-law fit of the form 

Var( qi ) = dL-^, (18) 
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FIG. 5: Log-log plot of [(q)], 
at several temperatures. 



as a function of system size L FIG. 7: Same as Fig. |^ but at temperature T - 



0.42. 
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TABLE II: Fits for Var(q ; ). Fit parameters for the fit in 
Eq. (^) for different temperatures. We cannot quote fitting 
probabilities since we have the same number of data points 
as variables. 



T 


a 


b 


c 


0.200 


0.00100 


0.0205 


2.55 


0.247 


0.00087 


0.0328 


2.76 


0.305 


0.00073 


0.0611 


3.16 


0.420 


0.00036 


0.1044 


3.40 



V. CONCLUSIONS 

To conclude, we have studied the low-temperature 
properties of the four-dimensional XY spin glass at low 
temperatures. Our main result is that the order param- 
eter distribution P(q) has, in addition to a peak, a tail 
that seems to extend for smaller values of q, and whose 
height seems to persist, as the system size increases, see 
Figs. U and |j. This interpretation of the data is compat- 
ible with the RSB picture or the KMPY scenario. How- 



FIG. 6: The distribution of the link overlap at T 
different sizes. 



0.20 for 



of the fits is poor as 
Q. The effective ex- 



see Table HI. However, the quali 
shown by the fitting probabilities!!, 
ponent fii is found to vary with temperature. Extrapolat- 
ing to T = 0, we obtain ^ = 6'+2(d-d a ) = 0.294±0.073 
If we assume that 8' w 0, this gives d—d s — 0.147±0.036 



TABLE III: Fit parameters for the fit in Eq. (jl|) for different 
temperatures. Note that the fit probabilities Q axe small. 
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Mi 
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0.200 


-4.92 ± 0.06 


1.07 ±0.05 


5.0 x 10" 2 


0.247 


-4.50 ± 0.05 


1.38 ±0.04 


3.6 x 10 -3 
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-3.99 ± 0.04 


1.77 ±0.03 


2.9 x 10 -6 
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fully equilibrated at low temperatures, but is rather only 
equilibrated up to some finite length scale, which only in- 
creases slowly with increasing measurement time. Thus 
a complete understanding will require a nonequilibrium 
theory, but a component of this is likely to be a theory of 
equilibrium on finite scales where local equilibrium has 
been achieved. 

We have also studied the link overlap qi . The variance 
of qi decreases with increasing L but we are unable to 
ascertain whether it tends to zero in the thermodynamic 
limit, and hence we are unable to determine whether or 
not the surface is space filling. 

In future work, it would be useful to look more care- 
fully at the nature of the large-scale low-energy excita- 
tions to see whether they correspond to gradual orien- 
tations in the spin directions or whether vortices play a 
role. 
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FIG. 8: Log-log plot of the variance of qi as a function of 
system size L at several temperatures. 
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